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Abstract 

Consider a large uniformly mixing dynamic population, which has constant birth 
rate and exponentially distributed lifetimes, with mean population size n. A Markovian 
SIR (susceptible —)• infective —)• recovered) infectious disease, having importation of 
infectives, taking place in this population is analysed. The main situation treated is 
where n —)• oo, keeping the basic reproduction number Rq as well as the importation 
rate of infectives fixed, but assuming that the quotient of the average infectious period 
and the average lifetime tends to 0 faster than 1/logn. It is shown that, as n —>■ oo, 
the behaviour of the 3-dimensional process describing the evolution of the fraction of 
the population that are susceptible, infective and recovered, is encapsulated in a 1- 
dimensional regenerative process S = {S{t)',t > 0} describing the limiting fraction of 
the population that are susceptible. The process S grows deterministically, except at 
one random time point per regenerative cycle, where it jumps down by a size that is 
completely determined by the waiting time since the previous jump. Properties of the 
process S, including the jump size and stationary distributions, are determined. 


1 Introduction 

The mathematical theory for the spread of infectious diseases has a long history and is by 
now quite rich (e.g., [Diekmann et al. (2013)| ). One of the more common type of disease 
models is called SIR (susceptible —>■ infective —)• recovered) meaning that individuals are at 
hrst Susceptible. If infected (by someone) they immediately become Infectious (being able 
to spread the disease onwards). After some time an infectious individual Recovers, which 
also means that the individual is immune to further infection from the disease. Such models 
were originally studied for populations assuming homogeneous mixing, but during the last 
few decades considerable effort has been put into analysing epidemic models in communities 
which are not homogeneously mixing but instead may be described using some type of social 
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structure, such as a community of households (e.g. [Ball et al. (1997)] ) or a random network 
describing possible contacts (e.g. |Newman(2002)| ). The vast majority of papers devoted to 
these type of problems assume a hxed community and community structure. 

In the current paper we treat the situation where the population is dynamic in the 
sense that people die and new individuals are born, or more precisely immigrate into the 
population. Further, we assume that there is also importation of infectious individuals 
(randomly in time according to a homogeneous Poisson process), implying that the disease 
never vanishes forever. In order to facilitate analytical progress we consider only the case of 
a homogeneously mixing community, which in network terminology corresponds to treating 
the complete network. 

Models for recurrent epidemics go back to the deterministic formulations of |Hamer(1906)| _ 

and Soper (1929)| . A stochastic treatment was given hrst in the pioneering work of |Bartlett(1956)| , 
who considered an SIR model with importation of both susceptibles and infectives, but with¬ 
out disease-unrelated deaths. An alternative model, with disease-unrelated deaths but no 
importation of infectives, has been studied extensively (e.g. |Nasell(1999)| and the references 
therein). Interest often centres on the time to extinction of infection and the closely-related 
problem of the critical community size for an infection to persist in a population. 

We consider a Markovian SIR epidemic with demography and importation of infectives, 
in which infectious individuals infect new individuals at constant rate and the infectious 
period is exponentially distributed. We study limit properties of the epidemic when the 
average population size n tends to inhnity. Our focus lies on the case where the limit is 
taken keeping the basic reproduction number Rq (i.e. the average number of susceptibles 
infected by a single infective in an otherwise fully susceptible population of size n) and the 
immigration rate of infectives hxed, whereas the quotient of the average infectious period 
and the average lifetime tends to 0 faster than 1/logn. For many infectious diseases this 
quotient typically lies between 10“"^ and 10“^, hence supporting this asymptotic regime, but 
in the discussion we treat other asymptotic regimes briehy. 

Under the above asymptotic regime, all epidemic outbreaks are short, having duration 
that tends to 0 in probability as n —>■ oo. Further, as n —>■ oo, epidemic outbreaks are either 
minor, having size of order Op{n), or major, having size of exact order 0p(n). It follows that, 
as n —)■ oo, the behaviour of the three-dimensional process describing the evolution of the 
fraction of the population that are susceptible, infective and recovered, is encapsulated in 
a one-dimensional regenerative process S = {S{t)]t > 0}, describing the limiting fraction 
of the population that are susceptible. During each cycle, the process S makes one down 
jump, corresponding to the occurrence of a major outbreak, and except for this increases 
deterministically, as minor outbreaks have no effect on 5^^ in the limit as n —)■ oo. (Here, 

5(A = t > 0}, where, for f > 0, with being the 

number of susceptible individuals in the population at time t.) Note that does not 
converge weakly to S in the Skorohod topology since the sample paths of S are almost 
surely discontinuous but those of almost surely contain only jumps of size n~^, so are 
close to being continuous. Thus to obtain rigorous convergence results, we consider two 
processes, and S^\ which coincide with except during major outbreaks during 
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which they sandwich S^"'\ and prove that both and 


Skorohod topology (Theorem 2.1). 


converge weakly to S in the 
It then follows that certain fnnctionals of converge 


weakly to corresponding fnnctionals of S (Corollary 

The paper is strnctnred as follows. In Section [2, we dehne the model and the limiting 
regenerative process, give an intnitive explanation of why S approximates for large 
n and present the main convergence results. In Section we derive some properties of 
the limiting regenerative process: the jump size distribution, the associated renewal time 
distribution and the stationary distribution. In Section]^ we present simulations supporting 
the convergence result and illustrating various features of the limiting process. In Section 
we prove the main results. We end in Section with a Discussion summarising our results 
and also exploring briefly additional questions, such as other asymptotic regimes. 


2 The epidemic model and main results 

2.1 The Markovian SIR epidemic with demography and importa¬ 
tion of infectives 

We now define the Markovian SIR epidemic with demography and importation of infectives 
(SIR-D-I). We consider the process to be indexed by a target population size n, which 
we assume is a strictly positive constant. The population model is an immigration-death 
process with constant immigration rate and linear death rate. For f > 0, let denote 

the population size at time t. Then increases at constant rate fin and decreases at 

rate /iA^^”^(t). The population size hence fluctuates around n, which is assumed to be large. 

The Markovian SIR-epidemic on this population is defined as follows. For f > 0, let 
and denote the number of susceptibles, infectives and recovered, re¬ 
spectively, at time t, so We assume that /•^"^(O) = 0 

and that —)■ sq as n —)■ oo, where sq G (0,1] is constant. (The value of i?^"'^(0) has no 

effect on the ensuing epidemic.) . A fraction of all births (i.e. immigrants) are infectives 
and the remaining births are all susceptibles, so births of infectives occur at rate fiUKn and 
births of susceptibles occur at rate /rn(l — k„). While infectious, any given infective infects 
any given susceptible at rate independently between each distinct pair of individuals. 

Thus, approximately, each infective makes infectious contacts at the points of a homogeneous 
Poisson process having rate A„, with contacts being with individuals chosen independently 
and uniformly from the whole population; a contact with a susceptible individual results in 
that individual becoming infected, while a contact with an infectious or removed individual 
has no effect. Each infectious individual recovers and becomes immune at rate implying 
that the infectious period is exponentially distributed with rate parameter y^. 

More formally, the process | /("'^(f), f > O} is a continuous-time Markov 
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chain, with state space and transition intensities given by 

(n) _ 

(n) _ 

(n) —1 \ 

1,1+1,r) “ ’ 

(n) _ 

*?(s,i,r),(s,j—l,r+l) 'yn'^i 

corresponding to birth of a susceptible, birth of an infective, death of a susceptible, death 
of an infective, death of a recovered, infection of a susceptible and recovery of an infective, 
respectively. 

We study specihcally the case where the average population size n tends to inhnity in 
such a way that 

(a) the total importation rate firiKn of infectives tends to a strictly positive constant fiK, 
so KnU —)■ K as n —)■ cx); and 

(b) the infection and recovery rates satisfy \nhn —)■ > 1 and 

A„/ logn —)■ cx) as n —)■ cx). 

For ease of exposition, we assume that n is an integer, so sequences of epidemic processes are 
indexed by the natural numbers. However, all of the results of the paper are easily generalised 
to the case of a family of epidemic processes indexed by the positive real numbers. 

To conclude, the parameters of the model are: n, the average population size; /r, where 
1/p is the average lifetime and fin is the population birth rate; Xn, the infection rate; 
where l/y^ is the average length of the infectious period; and Kn, the fraction of births which 
are infectious, so ^riKn is the birth (or importation) rate of infectives. 

2.2 The limiting process S 

Let t > O}, where is the “fraction” of the population 

that is susceptible at time t. The process S = {S{t);t > 0} can be viewed as the limit 
of as n —)■ cxD under the above asymptotic regime. It is a Markovian regenerative 
process (e.g. |Asmussen(1987)| , Chapter V), with renewals occurring whenever S{t) = I/Rq- 
Between each renewal S{t) increases deterministically according to the differential equation 

^'(t) = p(l-^(t)), (2.1) 

except for one down jump (from above 1/Ro to below l/i?o)- This implies that 

^(u) = 1 - (1 - l/i?o)e"^“ (2.2) 
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before the jump (if u denotes the time from the last renewal). The random time T from a 
renewal to the jump has distribution specihed by 


P(T < t) = 1 - exp 




1 - 


RoSiu] 


du 


{t > 0), 


with S{u) given by (2.2), so 


P(T < t) = 1 - -Ro + l) (t > 0). 


(2.3) 


The size of the jump is specihed by the value S(T—) of the process just prior to the jump. 
More precisely, S{T) = S{T—){1 — t{S {T—))), where for s > t{s) is the unique strictly 
positive solution to the equation (cf. [Diekmann et al. (2013)|, equation (3.15)) 


J- — q-Rqst 


(2.4) 


In epidemic theory r(s) is known as the relative fraction infected among the initially suscep¬ 
tible of an SIR epidemic outbreak in which a fraction s are initially susceptible and the rest 
immune. Hence, the size of the down jump is S{T—)t{S{T—)). After the down jump, S{t) 


increases deterministically according to the same differential equation (2.1) until the next 
renewal point, so 

S{T + t) = l-il- 5(T))e-^^ 0 < t < /i-i log[(l - 5(T))/(1 - 1/Ro)] 

and the inter-renewal time is T -|- log[(l — S'(T))/(1 — l/i?o)]- Illustrations of S are given 
in Section m 


2.3 Main results and heuristics 

We hrst explain heuristically why S can be viewed as the limit of as n —>■ oo under that 


asymptotic regime described in Section 2.1 Suppose that n is large. Then when no infective 
is present, all that happens is that individuals die and new ones are born at approximately 
the same rate /rn. Recovered (immune) individuals that die are replaced by susceptible 
individuals, so the fraction of susceptibles increases at rate /i(l — 5^"'^(t)) which explains the 
deterministic growth rate of S. 

After an exponentially distributed holding time, with rate parameter unKn ~ an 
infective is born into the community. If the fraction susceptible is below I/Rq, then 

the effective reproduction number = RqS^"''^ (t) is strictly less than one, implying that, with 
probability tending to one as n —>■ oo, a large outbreak will not occur, so continues 

to grow approximately deterministically. If (t) > 1 /Rq when a new born infective enters 
the community, then with approximate probability 1 — 1 /(RoS^'^^t)) that infective gives rise 
to a major outbreak that infects order Q(n) susceptibles (cf. [Diekmann et al. (2013)| , pages 
53 and 376), otherwise only a minor outbreak, which infects order o{n) susceptibles, occurs 
and (t) continues to grow approximately deterministically. This explains the distribution 
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for T, the time from a renewal until a down jump in S, which has time varying intensity 
given by /ik multiplied by the limiting major outbreak probability (cf. |Bartlett(1956)| ). 

If a major outbreak takes place, the size of the outbreak among the susceptibles is given 
approximately by t{S(T—))S^'^\T—) where S{T—) denotes the limiting (as n —)■ cx)) fraction 
susceptible just prior to the outbreak and r(s) is dehned above (cf. [Diekmann et al. (2013) , 
page 60). The duration of such a major outbreak is of order ©(logn/A^) (cf. |Barbour(1997) ) 
which tends to 0 by assumption. Thus, if there is a major outbreak it happens momentarily 
and, in the limit as n —)■ oo, the fraction susceptible after the outbreak, S{T), satishes 
S{T) = S{T-){1-t{S{T-)). 

Although the above heuristic argument makes it plausible that the normalised susceptible 
process converges to the regenerative process S, there are two complicating factors in 
making the argument fully rigorous. First, as explained in Section it is not true that 
5 (^ 1 ) S' as n —)■ 00 , where ^ denotes weak convergence in the space Zi)[0,cxD) of right- 
continuous functions / : [ 0 , cxd) — )■ M having limits from the left (i.e. cadlag functions), 
endowed with the Skorohod metric (e.g. Ethier and Kurtz(1986)| , Chapter 3). As explained 
also in Section 1, we overcome this problem by considering two processes, Si”^ and which 
coincide with except during major outbreaks, w hen they sandwich S^"\ and show that 
gM ^ ^ ^ ag ^ QQ- ggg Theorem 0 The second complicating factor is 

that the results referred to above concerning the probability, size and duration of a major 
outbreak are for an epidemic in a static population, whereas our population is dynamic. 
The results carry over to our setting because, in the limit as n —)■ 00 , the time scale of an 
epidemic outbreak is inhnitely faster than that of demographic change, but proofs need to 
be adapted accordingly. 

Before stating our main theorem, some more notation is required. Recall that 
is the number of infectives at time t in the SIR-D-I epidemic with average population size 
n and that we consider epidemics with no infective at time 0, i.e. with J("'^(0) = 0. Let 


4”^ = u'q’ = 0. For k = 1,2,---, let = inf{t > 




(n) 


U 


(n) 

k-1 


/(’^)(t) > logn} and 


(n) 

Ul = 


inf{t > = 0}. Thus, provided n is sufficiently large, the fcth major outbreak 


starts at approximately time and ends at time (The choice of logn to delineate 

major outbreaks is essentially arbitrary. Our proofs work equally well if logn is replaced by 
any function g{n) which satishes g{n) —)■ 00 and n~ 2 g{n) —)■ 0 as n —)■ cxd.) For t > 0, let 


(n) 


Ri")(t) = 


and 






lit ^ for some i. 


mm (n)< („) if t e = 1 , 2 , • • • , 




if t ^ for some z. 


max {„)< („) if t e nj”^), z = 1 , 2 , ■ ■ ■ . 

The following theorem is proved in Section |5.1 


Theorem 2.1 Suppose that lim„^oo 5'^"'^(0) = Sq. Then, as n ^ 00 , 

S and 4"^ ^ S, 
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where S'(O) = Sq. 


An immediate consequence of Theorem 2.1 is that suitable functionals of converge 


weakly to corresponding functionals of S. For g,h E Z1[0, cxd), let g < h denote g{t) < h{t) 
for all t > 0. A functional H : Z1[0, oo) —)■ M is called monotone if either Hf < Hg 
for all f,g E Z1[0, cxd) satisfying f < g, or Hf < Hg for all f,g E Z1[0, cxd) satisfying 
g ^ f ■ The following corollary, which can clearly be generalised to suitable non-real- 
valued functionals, follows immediately from Theorem 2.1 by using the continuous map¬ 
ping theorem (e.g. |Billingsley(1968)| ). For H : D[0,oo) —>■ M , let Ch = {f E Z1[0,cx)) : 
H is continuous at /}. 


Corollary 2.1 Suppose that lim^^oo = sq, H : Z1[0, oo) —)■ M zs monotone and 

P(^ e Ch) = 1. Then 

HsC) ffs as n —)■ cx), 

where S'(O) = sq. 


One functional which satishes the conditions of Corollary |2.1 is the hrst passage time 
functional Ha, dehned for given a E (0,1) by 

H f = S 

“ \ inf{f > 0 : f{t) < a} if /(O) > a. 

The functional Ha is clearly monotone and P (S' G Cna) = 1, cf. |Pollard(1984)| , page 124. 

Another functional satisfying the conditions of Corollary 2T is the occupancy time func¬ 
tional Hf,, dehned for any given C > 0 and a E (0,1) by 


HfJ = 


4/W<«} 


dh 


(2.5) 


This functional is again clearly monotone. The proof that P |^S G Ch^^J = 1 is given at the 
end of Section 15.11 


3 Properties of the limiting process S' 


We now outline some properties of the regenerative process S which can be obtained from 
renewal and regenerative process theory (e.g. [Asmussen(1987)], Chapters IV and V). As 


described in Section |2^ the stochastic part of the regenerative process is completely specihed 
by the waiting time T until the down jump, but it can be specihed equivalently by the jump 
size X = S{T-) - S{T). Noting that t{S{T-)) = (5(T-) - S{T))/S{T-)), it follows 


from (2.4) that 


S{T) 

S{T-) 


^ ^-Ra{S{T-)-S{T)) ^ ^-RoX 


7 



















whence 


S{T-) = 


X 


XqRoX 


and S{T) = 


X 


_ (3—RoX qRqX _ ^ j qRqX _ 2 ’ 

which can be used to obtain the distribution of the jump size X. The jump size is strictly 
less than r(l), as S{t) < 1 for all f > 0. Hence, for 0 < a; < r(l), 

Fx{x) = F{X < x) 


P fs(T- 


< 


X 


= P ( T < —jji ^ log 


^—Rox 


1 — x/{l — Q 


1 /i?o 


(using (2.2)) 


= 1 - 


- 1 

o 

1 

H 

1 

1 

o 

_1 

K 

O 

1 

(Ro-l)(l-e-^o^) 


1 — a; — e 


Rq 


(3.1) 


The lifetime distribution for the renewal process describing successive visits of S' to l/i?o 
may be derived as follows. During a cycle, the regenerative process S starts at 1/Rq and 


grows deterministically, according to (2.1), until the time T of the down jump. After this 


down jump it again grows deterministically, according to (2.1), until it reaches 1/Rq, when 


the next renewal occurs. If we change the order of these two parts, the process starts at 
S(T) and grows deterministically until it reaches S(T—). The lifetime T* hence equals the 


time it takes for the deterministic curve dehned by (2.1) to travel from SiT) to SiT—) This 
time equals 


T*=^l-Hog 


1 - S{T) 


= /i Mog 


^RoX 


1-X 


1-S{T-)J ^ ^ V(1-X)e«o^ - 1, 

This is a monotonic increasing function of X, so the renewal time distribution can be obtained 


numerically using the expression Fx{x) given by (3.1). 


The stationary distribution of S can be obtained using regenerative process theory 
(e.g. I Asmussen(1987) | , Chapter V, Section 3). During a regenerative cycle, the process 
S traverses s if and only if s lies between S{T) and S{T—). If it does, the density for the 
time spent there is inversely proportional to the derivative /i(l — s). Consequently, if we let 
fs* (s) denote the density of the stationary distribution of S, we have 


fs-(s) = 


ll(l - s) 


P(se [S(r),S(r-)]) (i - t(1) < s < i), 


(3.2) 


where c (= 1/E[T*]) is the normalizing constant making this a pdf. If s G [l/i^oA); then 
s G [*S'(T), S'(T—)] if and only if T > log((l — ^)/(l ~ ■^)- If s £ (1 l/-Ro); then 

s G [*S'(T), S'(T—)] if and only if X > g~^{s), where g : (0, r(l)) —)■ (1 — r(l), I/Rq) is dehned 


by g{x) = x/{e^°^ — 1). It then follows using (2.3) and (3.1) that, with s = g ^(s) 


fs*{s) = 


Ro(l-s-e-«o«) 

K 

(Ro-l)s 

(Ro-l)(l-e-«06) 


l—s—e~^0s 


Rfj 


1-6 
Rq —1 


kI 1 — 


Rq 


RqS^o 


if 1 — r(l) < s < 1/Ro, 
if 1/Ro < s < 1. 


































In the next section the density fs*{s) is calculated numerically and shown to agree with 
corresponding empirical values from simulations. 


4 Numerical illustrations 


We now present briefly some numerical and simulation results, which illustrate convergence 
of the epidemic process as well as properties of the limiting stationary distribution of the 
fraction susceptible S*. In Figure [T] the epidemic is simulated for 100 years in a population 
of n = 10, 000 individuals. In all hgures, Rq = 2 implying that the effective reproduction 
number = RoS^'^^t) is supercritical as soon as the population fraction susceptible exceeds 
I/Rq = 0.5. The average lifetime is l/fj, = 75 years and 7 = 50, so the average length of 
the infectious period is about 1 week. In the left panels of Figure [T| k = 20 , so the rate 
at which new infectives enter the population (/uk) equals 1 per 3.75 years, and in the right 
panels n = 200, so new infectives enter the population at rate 2| per year. The upper panels 
show the fraction of the population that is susceptible over the 100 period and the lower 
panels show the corresponding fraction that is infective. Observe that when k = 20 major 
outbreaks are less frequent but larger than when k = 200, and that there are appreciably 
more minor outbreaks when k = 100. Note also that epidemics are rarer than the importation 
rate of infectives suggests, for two reasons. First, major outbreaks can occur only when 
> l/i?o = 0.5, and secondly, when is above this threshold, major outbreaks 

do not occur each time an infective enters the community. In the lower left panel of Figure 
some minor outbreaks caused by importation of infectives can also be seen. 

In Figure realisations of the corresponding limiting processes are plotted. The same 
parameter values are used in both hgures. The stochastic features of the epidemic and the 
limiting process are in agreement, suggesting that the limiting behaviour has kicked in when 
n = 10,000. Note that, unlike in Figure there are no near-vertical lines as outbreaks are 
now instantaneous. 

We now illustrate properties of the stationary distribution of the fraction susceptible S*, 
both for the epidemic with n = 1, 000 and n = 10, 000, as well as for the limiting process. For 
the three processes, and for three different values of k, we simulate the epidemic and limiting 
processes for 10,000 years and in Figure]^ we plot bar charts of the relative time spent with 
specihed fraction susceptible. The processes are simulated over a very long time span so that 
the empirical distribution of the fraction susceptible is close to the corresponding stationary 


distribution. (Recalling the functional dehned at the end of Section 2.3, note that by 
standard regenerative process theory, for any hxed a G (0,1), —4 P(S'* < a) as 

as n —)■ cx).) The values of /i,7 and 
(Note that the value of 7, and hence also A (= i?o7)) is the 


f* —)■ CX) and, by Corollary 


2.1 


Rq are the same as in Figure [T 
same for both values of n.) The chosen values of k are k = 1, 3 and 100, corresponding to 
importation of infectious individuals on average one every 75, 25 and 0.75 years, respectively. 
In the plots we have also computed /s*(s), the stationary distribution of the limiting process, 
numerically as described in Section 

It is seen that the bar charts from the epidemics resemble the limiting stationary distri- 
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K = 20 K = 200 




K = 20 K = 200 




Figure 1: Simulation of the SIR-D-I epidemic with n = 10,000 individuals, Rq = 2. In the 
left panels k = 20 and k = 200 in the right panels. The average life length is l//i = 75 
years and mean infectious period is I/7 ~ 1 week. The fraction of the population susceptible 
(upper panels) and infective (lower panels) is plotted over a 100 year period in both cases. 
The dashed line in the upper panels shows the critical fraction susceptible so that the effective 
reproduction number Re = 1. Note that the scales for the fraction of the population infective 
are different in the two lower panels; major outbreaks are appreciably larger in the left hgure. 


10 























K= 20 


1 


K= 200 



qI-1-1-1-1-1 

0 20 40 60 80 100 

t (years) 


0.9 

0.8 

0.7 

0.6 



0.4 - 
0.3 - 
0.2 - 
0.1 - 

qI-1-1-1-1-1 

0 20 40 60 80 100 

t (years) 


Figure 2: Simulation of the limiting process S for the same parameter values as in the 
epidemics in Figure 


bution /5*(s), except when n = 1,000 and n = 100. When n is small, few outbreaks take 
place, so even if the outbreaks are large, the population fraction of susceptibles is close to 
1 most of the time, which explains why the stationary distribution S* is concentrated at 
values close to 1. For moderate values of k, the stationary distribution has positive mass for 
nearly all s values between 1 — r(l) = 0.2032 (the fraction susceptible after a major outbreak 
starting with the entire population being susceptible) and 1. The stationary distribution is 
seen to be concentrated around 1/Ro when k is large, owing to the fact that a new major 
outbreak occurs quite soon after the population fraction of susceptibles exceeds I/Rq, with 
the effect that the size of major outbreaks is generally small. These observations imply that 
the stationary distribution is not stochastically decreasing (nor increasing) in k. 


5 Proofs 


5.1 Proof of Theorem 2.1 


Let (fl, R, P) be a probability space on which is dehned a homogeneous Poisson process rj 
on (0, oo) having rate /i/t and let 0 < ri < r2 < • • ■ denote the times of the points in rj. For 


n 


= 1, 2, ■ • •, let 7]^"''^ denote the point process with points at 0 < r 


(n) 


< < 


where 




' = ■:^Tk {k = 1, 2, • • ■). Let denote the epidemic process indexed by n. Then 
gives the points in time when infectives immigrate into the population in We construct 
E(^) (n = 1, 2, • • ■) and S by hrst conditioning on rj. 

The process S is constructed as follows. Recall the dehnition of r(s) at 2.4 Between 


the points of r/, S{t) increases deterministically according to the differential equation (2.1). 


For k = 1, 2, ■ ■ •, S has a down jump to ^(rfc—)[1 — r(S'(rfc—))] at time r*. with probability 
max(l — (i?o>S'(rfc—))“^, 0) (independently for successive k), otherwise S continues to grow 


according to (2.1). Thus, S can be described as follows. Let R < t 2 < ■ ■ ■ be the times of 
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Figure 3: Bar charts of the relative time spent with fraction s susceptible for the epidemic 
(with n = 1,000 and n = 10,000) as well as the limiting process. Also plotted is the 
stationary distribution of limiting process fs*{s). Parameter values are: average life length 
equals l//i = 75 years, Rq = 2, mean infectious period I/7 ~ 1 week and k = 1,3 and 100. 
Bar charts are based on simulation over 10,000 years. 
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the down jumps of S, so these form a subset of the points of rj. Let 

f{x, t) = 1 — {1 — x)e“^* (0 < a: < 1 , t > 0 ), 


so, for hxed x, the solution of (2.1) with S'(O) = x is f{x,t). Let to = 0 and suppose that 
So = 5(0) is given. Then, for fc = 0,1, • • •, 


S{t) = f{sk,t-tk) (4 < t < 4+i), 


(5.1) 


where, for fc = 1 , 2 , ■ ■ ■, the initial value Sk = 4(1 — '?‘( 4 )), with 4 = 5(4—) = f{sk-i,tk — 
tk-i). The precise dehnition of the construction of (n = 1, 2, ■ ■ ■) is not relevant at this 
stage. 


We prove Theorem 2.1 by hrst proving the corresponding result for processes conditioned 
on rj. 


Lemma 5.1 Suppose that lim„_,.oo 5*'”'^(0) = Sq- Then, for T-almost all rj, 

5 i”^|r 7^5 and 5+”^|p^5 as n —>■ oo. 


(5.2) 


In order to prove Lemma 5.1, we need some more notation and an extra lemma (Lemma 5.2 
below). Recall that, for n = 1, 2, • ■ •, we assume /’^”^(0) = 0, that t^^^ = = 0 and that, 

for fc = 1 , 2 , ■ • •, tjf^ = inf |f > : I^‘^'>{t) > lognj and = inf > t]f'^ : /("■^(f) = 0 


For n = 1, 2, ■ ■ •, let Sg""^ = 5*^”^ for fc = 1, 2, • • •, let ' 

S^"‘\t) and = max S^^\t) 




c, = mm 


Lemma 5.2 Suppose that lim^^go 5*'"'^(0) = Sg. Then the following hold for F-almost all p. 
(%) Fork = l,2,---, u^f'^lp-^tk, 4”^|p-^4, Sk, 

(n)\D , _(n) I D ~ 

4 Ih —t Sk and 4 \V —t Sk as n —)• oo. 


(ii) For k = 0,1, ■ ■ ■, 


sup 


- / ( 4 "’.« 


u 


(n) 


\p 0 as n —)■ cx). 


(5.3) 


(in) 4 —t cxD as fc —)■ cxD. 


Proof. See Section 5.2 


□ 
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Proof of Lemma 5.1 First note that since results concerning convergence in distribution 
in the Euclidean space carry over in all essential respects to convergence in distribution 
in (see |Bilhngsley(1968)] , page 19), the Skorohod representation theorem implies that 
we may assume that the convergence in Lemma 5^ holds almost surely. Let A E IF he the 
set a; G H such that (i) for fc = 1, 2 ■ ■ ■ , 

lim u'^\uj) = tk{oj), lim = tk{oj), hm = Sk{oj), 

n^oo n^oo n^oo 


(ii) for fc = 0,1, ■ 


lim 




sup 


(n) 




u 


(n) 




= 0 ; 


(5.4) 


and (iii) tk{oj) —)■ cx) as A; —>■ oo. Then P{A\r]) = 1 for P-almost all rj. 

For g,h E Z1[0, oo), d{g, h) denotes the distance between g and h in the Skorohod metric 
(see lEthier and Kurtz(1986)], Chapter 3.5). Let g satisfy P{A\g) = 1. We show that for all 
uj E A, 

lim d{S^\u), S{u)) = 0 and lim d{S^\u), S{co)) = 0. 

n^oo n^oo 


It then follows that, under the Skorohod metric, both S^'^\g and converge almost 

surely to S, which implies (5.2). 

By Proposition 5.3 on page 119 of Ethier and Kurtz(1986)| , to show that d{gn,g) —t 0 
as n —)■ cx) it is sufficient to show that for each T > 0, there exists a sequence (A„) of strictly 
increasing functions mapping [0, cxd) onto [0, cxd) so that 


qn). 


lim sup |A„(t)—t|=0 (5.5) 

n^oo 0<t<r 


and 


lim sup \g{Xn{t)) - gn{t)\ = 0. 

n-s-oo o<t<T 


(5.6) 


For ease of exposition we now suppress dependence on oj. Fix T > ti and let m = 
max{fc : tfc < T}, so 1 < m < oo.. Then, there exists 5 > 0 such that u'm < T + 5 for all 
sufficiently large n. For such n, let A^^ be the piecewise-linear function joining the points 
(0, 0), • • • , {T +5, T +6), with = tiort > T +5. Similarly, let A+^ be 

the piecewise-linear function joining the points (0, 0), ■ ■ ■ , {u'm Am)-, {T + -|- 5), 

with X^ff\t) = t ioi t > T + 5. (Note that X^ff\f) < X^ff\t) with strict inequality for 
t E (0,T -|- 5).) The functions A^^ and A+^ are strictly increasing and satisfy (5.5), since 


that 


tF' —t tk and u 


-E- tk as n —)■ oo (fc = 1, 2, • • ■ , m). Thus, to complete the proof we show 


hm sup 

rn-oo o<t<T 


S{X^f!\t)) -S)f^’{t) 


(A I 


= 0 


(5.7) 
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and 


lim sup 

n^oo o<t<T 




(n)/ 


= 0 . 


Considering (5.7) first, note that for k = 1,2, ■■■ ,m, since S 
and = 4”^ for all t e 


(n) 


(5.8) 


IS increasing on 


sup 




(n) 


< max 




-'k 


—> 0 as n —)■ oo, (5.9) 

since = ^k, -t = 5(4) and a4^(m 4^) -)■ 4 (as -)■ tk and a 4^ is 

continuous), so S ^A4^(M4^)j —)■ Sfc as 5 is right-continuous. 

Also, for fc = 1, 2, • • • , m — 1, 


sup 


5 ( A4^(t)) — S'l’it) < A{n, k) + B{n, k). 


=,(n), 


where 


A{n, k) = sup S - / (4”^^ - “44 

fn^ ^ . fnl 5 / 5 / 




and 


S(n, k) = sup f (sP,t- 

Now a 444 ^ [4,4+i) for t e [m4^4+i)> so using ( |5A] ), 

A{n,k)= sup / 4fe,A444- 4 ) -/(s4\^-“44 

/ (^Sfc, A44^) - 4 ) - / (sfc,f - 


4"’<i<44i 


< sup 




-|- sup 




/ (^Sfc, f - m44 “ 4 (s4\ ^ - “44 • 


A simple argument using the mean value theorem shows that, for x G [0,1] and t, t' > 0, 

\f{x,t) - f{x,t')\ < (1 -x)/i|f-4|. (5.11) 


sup 

4"'<*<44i 


X^^\t) - tk - {t - 

< sup 

A44t)-t 

+ 

+ (A 

tk-ul’ 







—)■ 0 as n —)■ 00 , 
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as satisfies (5.5) and —)■ as n —)■ cxd. It then follows nsing (5.11) that the hrst 

term on the right hand side of (|5.10|) tends to 0 as n —)■ cx). Also, for x,y E [0,1] and t > 0, 


- f{y,t) = {y-x)e 


so the second term on the right hand side of ( |5.10[ ) tends to 0 as n —>■ cx), since —)■ Sk as 

n —)■ cx. Thus, A{n, fc) —)■ 0 as n —)■ oo. 


Note that S^\t)— = for t E t^k+i), so (5.4) implies that B{n,k) also con¬ 

verges to 0 as n —>■ cx, whence 


sup 




(n) 


—0 as n —)■ X. 


(5.12) 


Combining (5.9) and (5.12) yields that. 


lim 


sup 






W/ 


= 0 . 


(5.13) 


A similar argument to the derivation of (5.12) yields 


lim sup 

n^OO Cti') 

0<t<u[ ^ 


= lim sup 

k"^<t<T 




W/ 


= 0 , 


which together with (5.13) yields (5.7), as required. 


The proof of (5.8) is similar to that of (5.7) and hence omitted. 


□ 


Proof of Theorem 2.1\ We prove the result for S^ff\ The proof for is identical. Recall 
that if Xn (n = 1, 2, ■ ■ ■) and X are random elements of /1[0, x) then X^ X as n —)■ x if 
and only if E [/(X„)] —)■ [/(X)] as n —)■ x for all bounded, uniformly continuous functions 
/ : I1[0, x) —)■ M (see, for example, [Ethier and Kurtz(1986)| , Chapter 3, Theorem 3.1). Let 
/ : I1[0, x) —)■ M be any such function. Then Lemma 5.1 implies that, for P-almost all y, 


lim E 

n^oo 


/(si”>)|i;J =E|/(S)M. 

Hence, by the dominated convergence theorem. 


lim E 

n^oo 






= lim En 


E 


/(Si">)|>) 
ns'plr, 


= E, 


lim E 

.71^00 


= E,|E[/(S)|,|| 

= E[/(S)]. 

This holds for all bounded, uniformly continuous f : D[0, x) —)• M, so 
as required. 


S' as n —)■ X, 
□ 
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We end this subsection by showing that the occupancy time functional iff*, defined 
at (2.5), satisfies P G = 1- Recall that ti < t 2 < ■ ■ ■ denote the jump times of S. 

Let vi = inf{t > 0 : S{t) = a} and, for fc = 2, 3, • • •, let Vk = inf{t > Vk-i : S(t) = a}. Let 
C G -F be the set of a; G 12 such that tk{oj) (and hence also Vk{oj)) tends to oo as A; —)■ oo. 


Then, by Lemma 5.2 (hi), P(C) = 1. We show that if gn G Z2[0, oo) (n = 1,2, ••■) and 
lim„^oo d{gn, S{u;)) = 0, then lim^^oo Ht*9n = for u G C, whence P G = 

1 . 

Suppose that u E C. Dropping the explict dependence of S on u, since lim,^_,.oo d{gn, S) = 
0, by Proposition 5.3 on page 119 of [Ethier and Kurtz(1986)| , there exists a sequence (A„) 
of strictly increasing functions mapping [0, cx)) onto [0, oo) such that 


lim sup |An(t) 

n-^oo o<t<t* 


■t|=0 and lim sup |S'(A„(t)) — ( 7 „(f)| = 0. (5-14) 

n->-oo 0<t<t* 


Now 


where 


and 


\H-g^ - H-S\ = 


Rn — 


B„ = 


^{gn(t)<a} l{5(t)<a} df 




< An + Bn 


dt 


I l{S'(A„(t))<a} l{S(;:)<a} 


df. 


Let D = [0,f*] n ({ti, ^ 2 , • • ■ } U {ui, ^ 2 , • • ■ })• Then D has Lebesgue measure zero and 
l{S(A„(t))<a} - l{5(p< a} 0 as n CX), for t G [0,t*] \ D, since lim n^oo (t) = t, by 

the first equation in 5.14, and S is continuous at such t. Thus lim„_,.oo Rn = 0 by the 
dominated convergence theorem. A similar argument, using in addition the second equation 
in 15.141 shows that lim„_,.oo An = 0. Thus, lim„_,.oo Ht*9n = H^-S, as required. 


5.2 Proof of Lemma 15.21 


We prove Lemma 5^ by splitting the SIR-D-I epidemic process into cycles, where 
now a cycle begins at the end of a major outbreak and finishes at the end of the following 
major outbreak. Thus a cycle consists of two stages: stage 1, during which the susceptible 
population grows approximately deterministically until there are at least logn infectives 
present; and stage 2, comprising the major outbreak caused by these logn infectives, during 
which the susceptible population crashes. 

Recall that, as n —)■ oo, the point process describing immigration times of infectives 
in converges almost surely to the point process g governing times when down jumps may 
occur in the limiting process S. Lemma |5.4| considers the initial stage 1 and shows, using 


birth-and-death processes that sandwich the process of infectives, that for P-almost all g, as 
n —?• oo, for successive importations of infectives until a major outbreak occurs, the prob¬ 
ability a given importation triggers a major outbreak converges to the probability that the 
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corresponding importation results in a down jump in the limiting process S. Consequently, 
the time until there are at least log n infectives in converges weakly to the time of the 
hrst down jump in S, since converges almost surely to r]. Further, application of the 
law of large numbers for density dependent population processes ( [Ethier and Kurtz(1986)| , 
Chapter 11 ) shows that up until the hrst down jump of S, the scaled process of susceptibles, 
g{n) _ converges weakly in the uniform metric to S, since minor epidemics infect order 

Op{n) individuals. 


Lemmas 5^ and 5^ concern the limiting size and duration of a typical major outbreak. 
Lemma 15.51 considers outbreaks in which the initial number of infectives is of exact order 
n, for which the above-mentioned law of large numbers is applicable. This is then used to 
prove Lemma |5.6[ which considers major outbreaks trigg ere d by logn infectives. Finally, 


Lemma 5.2 follows easily by induction using Lemmas 5.4 and 5.6 


since 


£;(«) ig Markov. 


The proof involves extensive use of birth-and-death processes that bound the process of 
infectives in the epidemic model (cf. |Whittle(1955)| ). We hrst give some notation concern¬ 
ing birth-and-death processes and then state a lemma, proved in Appendix concerning 
properties of sequences of such processes. 

Let ^ > 0} denote a linear birth-and-death process, with = 

k, birth rate a and death rate (3. For x > k, let Ta^i3^k{x) = inf{f > 0 : Za^p^kit) > a:}, where 
Ta,/ 3 ,k{x) = 00 if Za^i 3 ^k{t) < X for alH > 0. (Throughout the paper we adopt the convention 
that the hitting time of an event is inhnite if the event never occurs.) Let Ta,i3,k{0) = inf{f > 
0 : Za,is^k{t) = 0 } denote the duration of Z^^js^k- For f > 0 , let Ba^^^kit) denote the total 
number of births during ( 0 ,f] in Za^p^k, and let Ba^p^k{oo) denote the total progeny of Za^^^k, 
not including the k ancestors. Further, for x > 0, let fa,i 3 ,k{,x) = inf{f > 0 : Ba^p^kif) > x}. 

Lemma 5.3 Suppose that a-a = a(3n {n = 1, 2, • • ■), where a > is constant and logn/ (3n —t 
0 as n —)■ cx). 


(a) If a < 1, then 

(i) for all t > 0, 


lim P(r„„,/ 3 „,i( 0 ) >t) = 0 ; 


(n) lim^^oo P (T'a„,/3„,i(fogn) = cx)) = 1; and 
(Hi) for any c> 0, 

( 0 ) ^ 0 as oo. 

(b) If a > 1, then 

(l) lim^^oo P (ra„,/3„,l(log^) < T-a„,/3„,l(0)) = 1 “ 
lim„^ooP(ra„,/3„,i(0) < T-„„A,i(fogn)) = 1 ; 

(ii) min(ra„,^„,i(logn),r„„,^„,i(0)) ^ 0 as n ^ cx); 


(Hi) lim„ 


P FB, 


(min{ra„,^„,i(logn),r„„,/ 3 „,i( 0 )}) < = 1 ; and 
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(iv) for any c> 0, 

(c'H') y 0 as n y oc. 

Before proceeding some more notation is required. For k = 1, 2, ■ ■ ■, let Xk = l{S(?'j,)< 5 (rj.-)} 
be the indicator function of the event that the kt\i point in rj yields a down jump in S. For 
n = 1, 2, • • ■ and k = 1, 2, • • •, let = inf |t > > logn or = o| and 

(n) _ .. 

Lemma 5.4 Suppose that ^^"'^(0) —^ sq as n ^ oo. Then the following hold for T-almost 
all p. 

(i) For k = 1,2, ■ ■ ■, 

lim P (xk^'' = 1 o.'ki'd = 0 for all i < k\p] 

n^oo \ / 

= P(xfc = 1 and Xi = 0 for all i < k\p). 


(ii) For k = 1,2, ■ ■ ■, as n ^ oo, 


sup 


1*5'*' \t) /(so,t)| for all i<k}\h ^ 0 - 


(in) For k = 1,2, ■ ■ ■, as n ^ oo, 


(’'*) 1 ^ 1 

and xl"'=o for all i<k} ^ and xi=o for all i<k} 


Proof. For ease of presentation we suppress explicit conditioning on p in the proof. First 
note that P (S'(ri—) = i?o = 0, since ri is a realisation of a continuous random variable. 
Assume without loss of generality that there is no recovered individual at time f = 0 . 
For f > 0, let Slf'\t) be the number of susceptibles at time t under the assumption that 
the immigration rate for susceptibles is fin{l — Kn) and the immigration rate for infectives 
is 0, and let S^\t) = SQ^\t)/n. Then, for any f > 0, application of Theorem 11.2.1 
of Ethier and Kurtz(1986)| (using the more general definition of a density dependent family 
given by equation (11.1.13) of that book) yields that, for any e > 0, 


lim P 

n^oo 


( sup 

- /(so,m) 

\0<u<t 




1 . 


(5.15) 


Recall that denote the epidemic process with average population size n. Con¬ 
sider the epidemic initiated by the immigration of an infective at time r^'^ in and let 
For ease of exposition, translate the time axis of so that the origin 
corresponds to r^\ With this new time origin, |J*^”)(t) : t > O} can be approximated by a 
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linear birth-and-death process t > 0 | having death rate 7 n + /U and (random) time- 

dependent birth rate given by XnSQ^\t). This approximation ignores depletion in the number 
of susceptibles owing to infection, so | t > o| is an upper bound for : t > O}. 

Let Si = /(so, ri) and fix e G (0, Si). Note that, with the change of origin, 5 q"'^( 0) = s["\ 


Then, using (5.15), for any 6 G (0,1), there exists i = t(e, h) > 0 and no = no{e, h) such that 


sup 
, 0<t<i 


- Si 


< - ] >1 -for all n > uq. 

2 - 2 - 


(5.16) 


For n > no and 0 < t < f , with probability at least 1 — the process |j*^"^(t) : t > 0 
is bounded below and above by the birth-and-death processes and ^Q,+ (e),g^i, 

respectively, where d“(e) = A„(si — |),Q!+(e) = A„(si -|- |) and /9n = 7n + l^- Further, 
a~{e)/Pn = Ro{si — for all sufficiently large n, the birth-and-death process 


since lim 


n^oo 


^o"(e),/ 3 „,i is bounded below by the birth-and-death process where a~(e) = Ro{si — 

e)l3n- Similarly, for all sufficiently large n, the birth-and-death process is bounded 

above by the birth-and-death process where a+(e) = i?o(si + e)/dn- 

Suppose first that Roh < 1. Then for all e G (0, Cq), where eo = Rq^ — Si, the birth-and- 


death process ^ is subcritical, so by Lemma 5.3 (a) (i), for all t > 0, 


lim P f r. 




(5.17) 


Setting t = t shows that, for all sufficiently large n, with probability at least 1 — 5, {/("■)(t) : 
t > 0}, and hence also {R'^'>{t) : t > 0}, is bounded above by ■^Q+(e) throughout its entire 
lifetime. Thus, 

lim inf P {x^i^ = o) > lim inf P (z +, s „ Nt) < logn for alH > o) —5 

n—>-00 V / n^oo \ )iP"i J 

= 1-5, 

by Lemma [T3| (a) (ii). Hence, since 5 G (0,1) is arbitrary, 

lim P (x^r^ = o) = 1 = P (xi = 0). 


Let = inf > 0 : > logn or = O}. Then it follows using (5.17) that 


D(^) 


0 as n —)■ oo. 


Suppose instead that RoSi > 1. Fix e G (0,ei), where ei = si — R^ , and 5 G (0,1). 
Then, similar to above, there exists ti such that, for all sufficiently large n, with probability 
at least 1 — |/^"'^(t) : t > o| is bounded above and below by -^a“(e)/ 3„5 respec¬ 
tively, throughout the interval [0,ti]. For x > 0, let = inf > 0 : /(’^)(f) > x|, 

fW( 0 ) = inf > 0 : = o| and = inf > 0 : > logn or = o|. 
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Note that the birth-and-death processes and -Z^+p) i are both supercritical. 

Then, by Lemma 5.3|^b)(ii), for all sufficiently large n, the process :t>0? is 


bounded below and above by ■Z^-(e) ji^-, respectively, throughout the interval [ 0 , 

Using Lemma [5.3[ b)(i), it then follows that 

liminf P (^^^"^(logn) < 

= l-vrTv^-(5.18) 

^o(5i-e) ^ ^ 


and 


limsupP (h^^^logn) < f(’^)( 0 )) 

72^00 

< limsupP (r„+(,)^^^^i(logn) < 5 

r?.—^no '' z 


= 1 - 


Ro{si + c) 


- 5 . 


Letting both e and 5 converge down to 0 in (5.18) and (5.19) yields 

1 


lim P (r^”'^(logn) < r^"^(0)) = 1 — 


RqSi 

= P(Xi = l). 


Further, using Lemma [5^b)(ii), it follows that 


(5.19) 


(5.20) 


f)in) 


0 as n —>■ cxD. 


(5.21) 


Recall that |/*^”^(t) : t > o| is an upper bound for |J*^’^)(t) : t > O}. We now show that 
the probability that the two processes coincide over converges to one as n —?• oo. 

In |j("^)(t) : t > o| births occur at time-dependent rate XnSQ^\t), whilst in : t > O} 

infections occur at time-dependent rate XnS^^\t). Now SQ"'\t) > for all f > 0, almost 

surely, so the two processes can be coupled by using an independent sequence Ui,U 2 ,--- 
of independent and identically distributed random variables that are uniformly distributed 
on (0,1), with the Rh birth in |j^"^^(t) : t > o| (which occurs at time U say) yielding an 

infection in {/*^”^(t) : t > O} if and only if Ui < /Sl"'\ti). 

For n = 1, 2, • • ■ and t > 0, let be the total number of births in |j("'^(t) : t > o| 

during (0,t]. Recall that the probability that |j(")(t) : t > o| is sandwiched between the 
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supercritical birth-and-death processes and throughout converges 

to one as u —)■ oo. It then follows using Lemma [5.3[b) (iii) that 


lim P > nA = 0. 

n^oo V / 

Also, since 0 as n —)■ cx), it follows using ( 5.16| ) that, for any e > 0, 

lim P (sAit) >h-e for all t e [0, = 1. 

n^oo V / 


(5.22) 


(5.23) 


Suppose that 5*^”^(£)("■)) < ns and, forhxede G (0,si), (t) > si —|forallf G [0,.D("')]. 

Then, S'^”^(fj) > sAiti) — ns, for i = 1, 2, • • • , so if pA denotes the probability 

that the ith birth in |j("')(f) : f > o| yields an infection in {/("■)(f) : f > O}, then 


whence 


(n) _ ^ ^ _ 

B(") (n)G)) 


ns 


sP(.U 


> 1 - 


_ 2 
n 3 

Si - e’ 


n 

2=1 


Pi 


(n) 


> 


> 


_ 2 
n 3 

Si — ( 

_ 2 
n 3 

Si - e 

_ 1 
n 3 


b (") (hG)) 


n3 


> 1 - 

Si - e 

—)■ 1 as n —)■ CX). 


Thus, recalling (5.22) and (5.23), the probability that {/("^(t) : t > O} and |j("')(t) : t > 0 
coincide over [0, converges to one as n —>■ cx, which, together with ( |5.20 ), yields 

hmPfxS") = l)=l = P(xi = l), 


and, together with (5.21), yields 


£}("■) ^L). Q ag 72 —)■ oo. 

We have thus proved parts (i) and (iii) for k = 1. Note that, since P {xk = 0 for all /c = 1, 2, • 
0 , when reverting to the original time axis, the probability that the total number of individ¬ 
uals infected during [0,w["'^] in is less than n^ tends to one as n —)■ x, which combined 




with (5.15) proves part (ii) when k = 1. Parts (i), (ii) and (iii) for k > 1 follow easily by 
induction since the processes {(S'^”)(f), Z<^”^(f)) : t > O} (n = 1, 2, • • ■) and S are Markov. 

□ 
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Before proceeding we state some well-known facts abont the final ontcome of the de¬ 
terministic general epidemic (e.g. Andersson and Britton(2000a)| Chapter 1.4). For f > 0, 
let s{t) and i{t) denote respectively the density of susceptibles and infectives at time t, so 
{s{t),i{t)) are determined by the differential equations 


ds 

dt 


—Rosi, 


di 

dt 


Rosi — i, 


(5.24) 


with initial condition (s(0),i(0)) = (sqAo), where s(0) > 0 and i(0) > 0. Note that time 
is scaled so that the recovery rate is 1. Then s(t) decreases with t, hmt_j.oo i(t) = 0 and 
hmt_^oo s(t) = Soo{so,io), where Soo(so,io) is the unique solution in (0,1) of 

e _ e ^-Ro(so-\-io—Soo) 

Oqq - 


Note that Soo is continuous in (sq, io) and Soo(so, *o) —t Soo(so, 0) as io 0, where (recall ( |2.4[ )) 

Sqo (SQj 0) 


0 

So(l - r(so)) 


if RqSo < 1, 
if RqSo > 1. 


In the following two lemmas, there is no importation of infectives in | t > O}, 

though births of susceptibles still occur at rate /in(l —k„). For t > 0, let 


Lemma 5.5 Suppose that (S'*^"'i(0),(so,*o) as n ^ oo, where So > ^ and 
zo > 0. Let = inf{t > 0 : = 0}. Then, as n ^ oo, 

(i) Soo(soAo), 

(^^) ^ 0 . 


Proof. For n = 1,2, and t > 0, let S^^\f) = S^'^^t/'jn) and I^'^\t) = 

Let = {xW(t) :t >0}, where X(")(t) = {s^^\t), . The process X(") is a 

continuous-time Markov chain with transition intensities 


(n) _ 


(1 ^n)/^ 


Ir, 




(n) 


(n) _ 




= n 


fi s 
Inn 



n 7 „n 


s i 
n n 


corresponding to a birth of a susceptible, a death of a susceptible, an infection of a susceptible, 
and a recovery or death of an infective, respectively. 
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The transition intensities are written in the above form to indicate that the family of pro¬ 
cesses : n = 1, 2, ■ ■ ■ } is asymptotically density dependent, as dehned by |Pollett(1990)| . 
Let E be any compact snbset of [0, cxd)^. Recall that «:„—)■ 0, 7 „ —)■ cxd and ^ Rq as 

n —>■ cx). Hence, as n —>■ cx), each of snp^^j ^^g^ snp( 3 , ,^)gg xy and 

sup(a:,j,)eE converges to 0. It follows that the conditions of Theorem 3.1 in |Pollett(1990)| 
are satished, whence, for any e > 0 and any t > 0, 


lim P 


( snp 

-X(")(t) -x(t) 

\0<u<t 

n 


<e =1, 


(5.25) 


where x(f) = is the solntion of the deterministic general epidemic (5.24) having 

initial condition (s(0),z(0)) = (so,^o)- Write Soo for Soo(so;*o)- There exists Cq > 0 snch that 
Ro{soo + ^o) < 1, since otherwise limt_,.oo *(t) wonld be strictly positive. Given e G (0,eo), 
choose e' > 0 so that 

/ Ro{soo + ^o) 


e 

< -. 


1 ~ Roi^oo + ^o) 8 


(5.26) 


There exists G > 0 snch that i{ti) < e' and s(ti) G [soo, s^o + f)- Then (5.25) implies that 


lim P 

n^oo 



1 and 


lim P 

n^oo 






1 , 


so, reverting to the original time scale and letting = H/yn, 

lim P ('|R(”^(t„) - Soo| < = 1 and lim P (^J^"^(t„) < = 1. (5.27) 

n^oo \' ' Z / n^oo \ ^ / 


Observe that, whilst -l- t) < Soo + e, the process +1) : t > O} is bonnded 

above by the birth-and-death process p^'n] ’ where = (soo + e)A„ and Pn = 1n + h- 

Now an/Pn —t Ro{soc + e) as u ^ oo, so, for all snfficiently large n, Z- ^ 
bonnded above by , where = Ro(soo + eo)/^n- 

Recall that Ba^y^k and Ta,i3,k{0) denote the total nnmber of births in and the extinction 
time of Za^y^k, respectively. Then 


E [Ba„,l3„,l] — 


Ro(Soo + ^o) 

1 ~ Roi^oo + Co) 


and, recalling (5.26), application of the strong law of large nnmbers yields 


lim P ( —B 


n^oo \ Tl 




< T = 1 


Also, Lemma 5.3[a)(iii) implies that 


T 


a„,l3„, [f e'n] 


36'nl (0) ^ 0 as n -)■ oo. 


(5.28) 


(5.29) 
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Recall that |5'g"'^(t) : t > o| denotes the process that describes the number of suscep- 

tibles in the absence of any infectives and suppose that Sq^\o) = S^‘^\tn)- For t > 0, let 
and be the total number of births and deaths, respectively, during (0,t] in 

|s'g"'^(t) : t > o|. Using (5.15) and the fact that B^\t) has a Poisson distribution with 
mean there exists t = t{e) > 0 such that 


ne\ 


lim P < — ) = 1 and lim P < ^ ) = 1. (5.30) 

n^oo \ 4 / n^oo \ 4 / 

The processes { + t), + t)) : f > O} and 0 [-e'n] coupled so 

that + t) < \^e'n \ (^) whilst +1) < Sqo + The hrst equations in (5.27) 

and (5.30) imply that 


lim P ( sup + t) < Soo + e 1 =1, 

. 0<t<i 


so (5.29) implies that, with probability tending to 1 as n —)■ cx), the coupling holds thoughout 
the lifetime of Z^ ^ ‘ ^ inf{f > 0 : = 0}. The coupling implies 


that — tn < « rs /„i, so part (ii) of the lemma follows from (|5.29|), since —)■ 0 

(^ntPni 2^ ^ ^ ^ 1 I 

as n —)■ cx). Further, is at most the sum of S^"'\tn) and the number of births 

in and at least the difference between S^"'\tn) and the sum of the number of 

susceptible deaths in and B^^^^ [Se'n] > ( |5.27| ), (5.28) and (5.30) imply that 


lim P (^Soo - e < < Soo + = 1, 

n^oo y 4 / 


proving part (i) of the lemma, since e G (0, eg) can be arbitrarily small. 

Lemma 5.6 Suppose that /("■^(O) = [logn] (n = 1, 2, ■ ■ ■) and 

^("■)(0) -hy sg Qg yi ^ oo^ where Sg > ^. Let = inf {t > 0 : I^'^^t) = O}, 


□ 


(n) 

c\ = mm 


S^'^\t) and = 


max 

0<t<u["^ 


S^^\t). 


Then, as n ^ oo, 

(i) so(l — 'r(so))? where the function r(s) is defined at (2.4) 

(^^) 0 ; 

(Hi) ^ so(l — t(so)) and ^ sq. 
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Proof. Fix 6 G (0,1) such that So(l — 36)Rq > 1. Then, whilst > So(l — 26), 

: t > 0} is bounded below by the birth-and-death process ^a„( 6 i),/ 9 „,|'iogn], where 
^n{6) = AnSo(l — 29) and /3n = 7n + P- Now an{6)/(3n —t -RoSo(l — 26) as n —)■ cx), so, 
for all sufficiently large n, is in turn bounded below by Za„{e),p„,\\ogn\, where 

Otnif^) = -Ro'So(l ~ ^^)Pn- 

Recall that denotes the epidemic process indexed by n. For t > 0, let B^^\t) be 
the total number of infections in during (0,t]. Let = inf > 0 : B^'^\t) > 6*Son}. 
Dehne : t > o| as in the proof of Lemma 


5.4 


For f > 0, let B^\t) and Djf^\t) 


be the total number of births and deaths, respectively, during (0,t] in Sl^\t) : t > 0 [>. As 


at (5.30), but note that S'q"'^( 0) is different here, for any e > 0 there exists i{e) > 0 such that 
lim P (Blf'\i{e)) < = 1 and lim P = 1. (5.31) 

n^oo \ 4 / n^oo \ 4 / 

Also, since ^^”i(0) So as n —)■ oo, 

lim p(|^(’^)(0)-so| < = 1. (5.32) 

n^oo \ Z/ 

Observe that, if < £(e), |S'*^"'i(0) — So| < Bjf^\i{e)) < ^ and Dlf\i{e)) < then 


max . < So(l - 0) + ye, 




(5.33) 


obtained by making S'^"i(0) and Blf^\i{e)) as large as possible and assuming no susceptible 
dies during [0,rg"'^], and 


mm 


0<t<T, 


(n) 




(5.34) 


obtained by making 5i'‘i(0) as small as possible, Dg\i{ej) as large as possible and assuming 
no susceptible is born during [0, ]. 

Recall that, whilst > so(l — 29), {/("'^(t) : t > O} is bounded below by the birth- 

and-death process Za„{e),i 3 „,\iogn], so < ^ar!,i 3 „,\\ogn-]i^o 6 n), provided S^^\t) > so(l - 26) 
throughout [0,Tg^'^]. Now 13^ ^iogn]i^o^'^) 0 as n —>■ cxd, by Lemma 5.3'b)(iv), so 


P f 


'a„ pn [logn] ^ —t 1 as n —>■ OO, for any e > 0. Setting e = so 6 in (5.34), 


lO) 


using 


(5.31), ( |5.32 ) and noting that n sq 6 as n —)■ 00 , shows that 


lim P 

n—1-00 


mm 


, 0<t<T, 


(n) 


>so(l-20) =1, 


so 


r. 


(n) 

e 


0 as n —)■ 00 . 


(5.35) 
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Further, since for any e > 0, P < t(e)j —)■ 1 as n —)■ oo, it follows from (5.31 )-(5.34) 
that, for any e > 0, 


lim P 

n^oo 




so 


Sq{1 — 9) as n —>■ oo. (5.36) 

It is straightforward to couple the jump processes of : f > O} and 

/ \ / \ f ^ 

{■^Q!„(0),^n,|'logn] (^)) .Sa„(6»),/3„,|'logn] (^) • ^ ^ 0| tO shoW that I [Tq ) > ■Z^a„(e),/ 3 „,|'logn] (Fg ), 

/ \ / \ st 

where Tq = [logn] and > denotes stochastically greater than. Further, recalling 

that ^a„(e),/ 3 „,riogni h as the same distribution as {^Ro(i- 30 ),i,riogni(/dn^) : ^ > O}, it follows 
using |Nerman(1981)| , 

Theorem 5.4, that 


7 f'-M] 

^an(d),/3n,[logn'][Tg ) a.s.^ 


1 - 


-^ 050(1 ~ 30 ) 


as n oo. 


-Sa„(0),/3„,|'logn] IF? . 

Thus, since n-^Ba„( 0 )j„,[iogn](fg'"^) ^ So0 as n oo, 

lim P >i-(0)) = 1, 

n^oo \ J 

where 


(5.37) 


i.{d) = 


1 - 


Sq^^. 


-Ro'So(l ~ 36) 

A similar argument using an upper bounding birth-and-death process yields that 


lim P ' 

n^oo 


where 


hW = 


< 4W) = 1, 
2 


(5.38) 


1 + 


SqO. 


-Ro'So(l — 36) 

Exploiting the Markov property of {(S'^"'^(f),/*^”^(f)) : t > O}, ( |5.36 )-(5.38) and Lemma 


5.5 


imply that, for any e > 0, 

lim P (soo(so(l - 6),i_{6)) - e < < Soo(so(l - 6),i+{d)) + e) = 1. 

Letting 6 I 0, noting that z_(0+) = i+(0+) = 0 and using the continuity properties of Soo, 
yield that, for any e > 0, 


lim P 

n^oo 


5(”^(mS”^) - so(l - r(so)) 


<6 = 1 , 


proving part (i) of the lemma. Part (ii) follows immediately using (5.35), 

(5.36), (5.38) and Lemma 5A|(ii). Part (hi) is an easy concequence of parts (i) and (ii) 
and (5.31). □ 


27 


































P roof of Lemma 5.2 The l emm a follows easily by induction using Lemmas 5.4 and Lem¬ 


mas 


5.6 


Lemma 


First note Lemma 


5.4 


5.4 


(i) and (iii) imply that ti as n —>■ oo, and 

(n) 


D 


(i) a nd ( ii) imply that (5.3) holds for A; = 0 and —)|n 




n —>■ oo. Lemma |5.6| (ii) then yields tha 
that Si as n —)■ oo, and Lemma 


u 


(n) 


5.6 


\t] 


D 


as n —>■ oo, Lemma 


as n —)■ oo. Now {t > O} is Markov, so, since si as 

n —)■ oo, the above argument can be repeated for fc = 2, 3, • • ■. Part (iii) is immediate, since 
{Ai, ^ 2 , • ■ ■ } ^ {^ 1 ) ^ 2 ,- " where ri, r 2 , • ■ ■ are the times of the points in rj. □ 


5.6 


-)■ Si as 
i) yields 

iii) yields that si and si 

S(n)l 


6 Discussion 

In the paper it is proved that for an SIR epidemic in a dynamic population (whose size 
fluctuates around n), in which there is importation of infectives at a constant rate, the 
normalised process of susceptibles converges to a regenerative process S' as n —>■ oo. Further, 
properties of the limiting process S are derived. The asymptotic regime considered is for the 
situation when the rate of importation of infectives Kfa and the basic reproduction number Rq 
remain constant with n, whereas the average length of the infectious period l/ 7 n converges 
to 0 faster than 1 / log n (in most real-life epidemics, the ratio of average infectious period 
and average lifetime lies between 10 “^ and 10 “^). 

Other asymptotic regimes could of course also be considered. For example, if the im¬ 
portation rate of infectives grows with n, then there will always be infectives present in the 
population resembling an endemic situation. If the duration of an infectious period remains 
hxed (or at least grows slower than logn), then the duration of a single outbreak will be 
long and the typical time horizon will not go beyond the first outbreak. A more complicated 
and interesting scenario seems to be for the asymptotic situation treated in the current pa¬ 
per, but where the epidemic is initiated with a fraction 1/Rq of the population susceptible 
and a large enough number of infectives. It then seems as if an endemic equilibrium will 
stabilize, but determining and proving this rigorously remains an open problem. For large 
but hnite n, it is possible for the process to get stuck in an endemic situation near the end 
of a major outbreak (with states similar to those just described). Eventually the epidemic 
leaves this endemic state and returns to the behaviour of the limiting process. In Figure]^ 
such a simulation is presented. The parameter values are n = 100, 000, /i = 1/75, k = 1 (so 
the importation rate of infectives is one per 75 years), Rq = 2 and 7 = 2 (so the average 
infectious period is 6 months). The left and right plots show the fraction of the population 
that are susceptible and infective, respectively, as functions of time. A quasi-endemic phase 
lasts roughly from years 1, 300 to 3, 000. Observe that major outbreaks become smaller prior 
to the process entering the quasi-endemic phase and fluctuations in the number of infectives 
increase in amplitude prior to the end of the quasi-endemic phase. Beside studying other 
asymptotic regimes, it could be of interest to increase realism in the model, for example, 
by relaxing exponential distributions of infectious periods and lifetimes and allowing for a 
latent state (cf. [Andersson and Britton(2000b)|, who consider epidemics with importation 





















Figure 4: Plot of an epidemic exhibiting quasi-endemic behaviour. 


of susceptibles only) or by having some population structure, such as network or households 


(see the challenges in |Pellis et al. (2015)| and |Ball 

et al. (2015)1). 


A Proof of Lemma 5.3 




Suppose that a < 1. Then, for f > 0, 




lim P(r«„,^„,i(0) > f) = 

lim P 

n^oo 

IV 


< 

lim E [Za^ 

n^oo 



= 

exp (—(1 — 

a)/ 3 nt) — !• 0 as n —)■ oo. 

since /?„ —)■ oo as n —)■ oo, proving part 

(a)(i). 



Observe that, for any k, 

IV 

lib 

{Za,l,ki./ 3 nt) : t 

> 0 }, where = denotes 

equal in distribution. It follows that 

./3n,l(l0gn) = 

^ra,i,i(logn). 

so 

lim P(r„„,y 3 „,i(logn) = 

n^oo 

oo) = lim P (ra,i,i(logn) 

n^oo 

= oo) 


= 1 , 

since is subcritical, proving part (a)(ii). 

For any t > 0, 


P {'Tan,l3ri,\cn] ^ P (^) ^ l) 

< E (t)] 

= \cn\ exp (—(1 — a)(3nt) 
—>■ 0 as n —)■ oo, 

since logn//3n —)• 0 as n —)■ oo, proving part (a)(iii). 
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Suppose that a > 1. Note that, since {Za„,i 3 „,i{t) : t > 0} = {Za^i^i{/3nt) : t > 0}, coupled 

realisations of : n = 1,2, • • • } can be obtained by setting Za„,i 3 „,i{t) = Za^i^i{f3nt) 

{n = 1,2, > 0). Now, see e.g. [Athreya and Ney(1972)| , page 112, there exists a 

random variable hh > 0, satisfying W = Q ii and only if Za,i^\ goes extinct, such that 

e-(“-i)*Za^i^i(t) W as t —>■ oo. (A.l) 

Note that, for any x > 0, = /l“^ra^i^i(x). If hh = 0, then ra^i^i(O) < oo, so 

T'an,/3n,i(0) i 0 as n —)■ cx), and maXi>o i i(f) < oo, so ra^i_i(x) = oo for all sufficiently large 
X, whence i(logn) = oo for all sufficiently large n. If W > 0, then ra^i^i(O) = oo, so 

'ra„,/3„,i(0) = oo. Also, for any f > 0, 


Za^,p^,i{t) > logn 


ZaXlif^nt) > logn 

e-(“-i)/^"*Z<,,i,i(/l„f) > e-(“-^)^"*logn. 


Now logn —)■ 0 as n —)■ oo, since lim„_,.oo log?7,//l„ = 0, so, since 11^ > 0, (A.l) 

implies that Za„^p^,i{t) > logn for all sufficiently large n. This holds for any f > 0, so 
'ra„,/ 3 „,i(logn) —)■ 0 as n —)■ oo. Part (b)(i) follows since P(IT = 0) = „ (t')(ii) also 

follows, indeed we have shown that, under the coupling, min (rQ„^/ 3 „ 4 (logn), 0 

as n —)■ oo. 

Further, it follows using |Nerman(1981)|, Theorem 5.4, that 


^ as t ^ oo. 


a — 1 


(A.2) 


where IF is the same random variable as in (A.l). If IF = 0, then i?ai^i(oo) < oo, so 
^an,/ 3n,i(T 'a„,/3n,i(0)) = -Ba,i,i(oo) < U3 for all Sufficiently large n. If IF > 0, then ( |A.lD 
and (A.2) imply that liiut^oo -Ba,i,i(^)/-^a,i,i(^) = so, since ra^i^i(fogn) — >■ oo as n — )■ oo, 

hm,^^oo5a„,/3„,i(fo„,/3„,i(logu))/fogn = whence 5„„,;3„,i(r„^,^„,i(logn)) < for all 
sufficiently large n. Part (b)(iii) now follows. 

A similar argument to the above shows that, if IF >0, then for any c > 0, Tar,,(in,i{cn) —)■ 
0 as n —)■ oo. Thus TQ,^,/ 3 „,|'iogn] (cn) 0 as n —>■ oo, since -^a„,/ 3 „,|'iogn] is the sum of 

st 

[logn] independent copies of Part (b)(iv) follows since jiogn](cn + 1) > 

Xa„,/3„,riogn] (cu). 
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